arXiv:1506.06398vl [astro-ph.CO] 21 Jun 2015 


Draft version June 23, 2015 

Preprint typeset using D-T^;jX style emulateapj v. 03/07/07 


FAST MAGNETIC FIELD AMPLIFICATION IN THE EARLY UNIVERSE: GROWTH OF COLLISIONLESS 

PLASMA INSTABILITIES IN TURBULENT MEDIA 

D. FALCETA-GONgALVES^’^ & G. KOWAL^ 

Draft version June 23, 2015 

ABSTRACT 

In this work we report a numerical study of the cosmic magnetic field amplification due to colli¬ 
sionless plasma instabilities. The collisionless magnetohydrodynamic equations derived account for 
the pressure anisotropy that leads, in specific conditions, to the firehose and mirror instabilities. We 
study the time evolution of seed fields in turbulence under the influence of such instabilities. An ap¬ 
proximate analytical time evolution of magnetic field is provided. The numerical simulations and the 
analytical predictions are compared. We found that i) amplification of magnetic field was efficient in 
firehose unstable turbulent regimes, but not in the mirror unstable models, ii) the growth rate of the 
magnetic energy density is much faster than the turbulent dynamo, iii) the efficient amplification oc¬ 
curs at small scales. The analytical prediction for the correlation between the growth timescales with 
pressure anisotropy ratio is confirmed by the numerical simulations. These results reinforce the idea 
that pressure anisotropies - driven naturally in a turbulent collisionless medium, e.g. the intergalactic 
medium -, could efficiently amplify the magnetic field in the early Universe (post-recombination era), 
previous to the collapse of the first large-scale gravitational structures. This mechanism, though fast 
for the small scale fields (^kpc scales), is however unable to provide relatively strong magnetic fields 
at large scales. Other mechanisms that were not accounted here (e.g., collisional turbulence once 
instabilities are quenched, velocity shear, or gravitationally induced inflows of gas into galaxies and 
clusters) could operate afterwards to build up large scale coherent field structures in the long time 
evolution. 

Subject headings: intergalactic medium — magnetic fields — turbulence — MHD 


1. INTRODUCTION 

The baryonic fraction of the intergalactic medium 
(IGM) is composed by dillute warm/hot plasmas. De¬ 
pending on the local environment the plasma proper¬ 
ties of the IGM may be different. For instance, the 
most diffuse component of the IGM, known as inter¬ 
galactic “voids”, are vast regions between clusters of 
galaxies in the Universe as seen today. Filaments of 
somewhat compressed and warmer (10® — lO’^K) gas 
form what is broadly understood as the IGM, while the 
denser (^ 10“®cm“®) and shock heated (^ 10^ — 10®K) 
plasma in clusters of galaxies are known as the intra¬ 
cluster medium (IGM). For many decades the intraclus- 
ter medium is known to be magne t ized (see reviews by 
Kronber^ll994 lEnfilin et ^l2005l IWidrow et al.l [20T4 
Durrer fc Neronmll201.‘ll and references therein), but un¬ 
fortunately any contraints on the magnetization of the 
IGM and intergalactic voids lack of more substantial ob¬ 
servational confirmation. The magnetic intensity and its 
respective correlation length must naturally be related to 
the processes from which they originate and/or amplify, 
which depend on local dynamical and physical proper¬ 
ties. However, the origin and/or amplification of the 
magnetic fields at cosmological scales are not completely 
understood yet, and represent a long standing major is¬ 
sue in modern astrophysics. 

Seed fields may have been generated at the pre- 
equipartition epoch by subatomic scale processes, such as 
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hadron phase transitions, or at later (post-equipartition) 
and pre-galactic epoch by large-scale magnetohydrody- 
namical (MHD) process es, e.g. the Vp x Vp battery 
term from Ohm’s Law (jBiermannI 1195011 , currents in¬ 
duced by electronic scat t ering of anisotropic rad iation 
fields (|Lanefer et al.ll200^ iDurrive fc Langeill20I4ll , evo- 


lution of turbulent fields le.g. lBaneriee fc .Iedamzikll2004 


and ma ny others), and even shock excited Weibe l insta ¬ 
bilities (ISchlickeiser fc Shuklal l2003t ISchlickeise^ l2005f) . 
(|.Iedamzik et al.lll998l l suggests however that photon dif¬ 
fusion damping would affect post-equipartition, but pre¬ 
recombination, magnetic fluctuations. Therefore most of 
the amplification processes mentioned above would only 
be important after recombination era. Theese could in 
principle generate magnetic fields with intensities of the 
order of 10“®° — 10“^®G. There are several implications 
of a strong magnetic field on the cosmological evolution 
of the Universe, from primordial nucleosynthesis to the 
statistics of the cosmic microwave background (CMB). 
Gurrent data limits the primordial (comoving) magnetic 
field intensity to Hiiunr < 5nG at IMpc lengthscale (see 
iTrivedi et al.il20fol : iPianck Collab .12014 for references on 
WMAP and PLANCK data analysis, respectively). Such 
field intensities are belo w the values inferred for the IGM 
{B ^ 0.1 — 10 pG fsee lGovoni fc Ferettill2004 and ref¬ 
erences therei n)), as well as for more diffus e groups of 
galaxies (e.g. iNikiel-Wroczvnski et al.ll201,'ll) . Recently, 
arguable lower limits for the magnetization of the diffuse 
intergalactic medium have been also p rovided based on 
gamm a-ray observations of blazars fsee lNeronov fc VovkI 
l2nTnli . These evidences combined together suggest that 
the Universe was magnetized even before the structure 
formation. Therefore, the problem seems to be twofold: 
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the first being the actual origin of the magnetic fields 
(seed fields), which may have occured very early in the 
Universe history (pre-equipartition); and the second be¬ 
ing the amplification of such seed fields to the values 
observed in the diffuse gas of clusters of galaxies. 

Observations of the Faraday rotation (FR) effect pro¬ 
vide currently the best estimates for the IGM mag¬ 
netic fields. In the FR effect the position angle of 
linearly polarized radiation shifts on the plane of sky 
as a function of wavelength A as A(j) = RM x A^, 
with the rotation measure RM being i?M (rad/m^) ~ 

812 ne{cm~^)B\\{fj,G)dl, where |1 denotes the di¬ 

rection parallel to the line of sight. 

It is well-known that RM strongly depends on the spa¬ 
tial distribution of the electron density (rie) and the mag¬ 
netic field geometry along the line of sight (i3||), and basic 
geometrical assumptions with respect to them are usually 
made. It is also natural to presume that the IGM mag¬ 
netic field is not uniform. Faraday rotation maps present 
ordered magnetic fluctuations at scales from ^ Ikpc as¬ 
sociated with cooling flows up to ^ IMpc associated to 
powerful active galactic nuclei (A GN) jets and cluster 
radio relics le.er. iFerreti et al.ll20I2tl . These lengthscales, 
however, cannot be directly understood as correlation 
scales for the magnetic field. Observational estimates 
of typical lengthscales usually assume an ad hoc corre¬ 
lation between the magnetic field and local plasma den¬ 
sity in the range o f 1 — SOkpc (e.g. iBonafede et al.l[201Cll : 
IVacca et al.l 1201 2 tl . Unfortunately, given the turbulent 
nature of the IGM, the actual correlation length (/corr) 
of the IGM/ICM magnetic field depends also on the spa¬ 
tial distribution of the velocity field fe.g. iBu rkhart et al.l 
120091 : IXu et al.l[200^ iFalccta-Goncalves et al.ll2ni4l l. and 

4orr for the magnetic field of the IGM cannot be ob- 
servationally estimated yet. Still, even though around 
equipartition level with the thermal and kinetic counter¬ 
parts, the presence of strong magnetic fields {B > fiG) is 
surprising, given the general understanding that /xG scale 
fields would be the result of galactic dynamo amplifica¬ 
tion of much less intense pre-galactic seed fields and the 
related timescales. 

During the last decade several authors have employed 
numerical simulations and analytical approximate so¬ 
lutions of collisional plasmas to study the magnetic 
amplification due to turbulent dynamo, galactic fields 
diff used into the ICM by outflows, and AGNs (see 
e .g. iDolag et al.ll 2002 HDubois et ani2009l: IPonnert et ahl 


2009t FaIceta-Gonca.lves et al.l 120101: Xu et al.l ^121: 


a-tjoncaives et al.l . , , 

Schober et al.l l2012t 120131 : IChol 12014 iFederrath et ahl 
2014 and references therein). It was found that all these 
processes could in principle account for the amplification 
of a pre-existing seed field to observable values. Such de¬ 
generacy among different mechanisms could be removed 
if the correlation lengths of density, velocity and mag¬ 
netic fields were directly and independently obtained ob- 
servationally. This because the different processes sug¬ 
gested so far operate at either different scales, or with 
different cross-correlations between the physical param¬ 
eters (B,n,p,v). However, the possibility that mag¬ 
netic fields could be present in vast volumes of the more 
diffuse intergalactic medium, where the density of galax¬ 
ies and the impact of AGNs are relatively small, points to 
the turbulent dynamo as a viable and ubiquitous process. 


Great effort on the understanding of turbulent 
dynamos has been employed by means of analyt¬ 
ical a nd numerical stud ies _ fse e ISch ekochih in et ahl 

2007t [Brandenburg et al.l I20l4 ISchober et ahl 120121 : 
de Gouveia Dal Pino et al]l201.'ll and references therein). 
The turbulent dynamo is the process by which kinetic 
energy of turbulent motions is converted into magnetic 
energy. Turbulent cells stretch and fold field lines in 
a weakly magnetized plasma increasing the total mag¬ 
netic energy. Such a dynamo can in principle operate 
at the large range of scales where there is turbulence, 
depending on t he propert i es of the flow. This was first 
pointed out by iBatcheloil (|195r)tl , and f urther developed 
into the Kraichn an-Kazantsev theory (iKazantsevI 119671 : 
lKraichn^ll968ll . Here, and throughout this paper, we 
assume the intergalactic plasma to be a highly conduct¬ 
ing medium in the sense that the resistivity ( 77 ) is much 
smaller than the viscosity (v), i.e. the Prandtl number 
Pm = I'/r] ^ 1. In this case, the small-scale dynamo 
(SSD) theory predicts an exponential growth of mag¬ 
netic energy, up to equipartition at the viscous scale. 
The power spectrum of the magnetic fluctuations is then 
proportional to resulting in a concentration of en¬ 
ergy at small scales. Once saturated at small scales, the 
dynamo enters in a nonlinear phase showing slower am¬ 
plification rates (linear with time). The scale at which 
turbulence is dissipated regulates the efficiency of non¬ 
linear phase for the dynamo process. If we consider the 
standard collisionality of the IGM, the non-linear phase 
of the turbulent dynamo is possibly of little interest for 
cosmological magnetic field studies. It has been sug¬ 
gested though that the enhanced effective collisionality 
due to plasma instabilities could even completely sup¬ 
press the kinetic phase, and the nonlinear dynamo woul d 
be all that is left (see iMogavero fc Schekochihii] I20l4l . 
As we discuss below, for the purpose of this work, the 
exact determination of the transition between kinetic^ 
and nonlinear regimes is not relevant. Due to the little 
understanding of the anomalous collisionality in diffuse 
media, and the fact that the nonlinear regime is even 
slower than the kinetic one, we will concentrate only on 
the early kinetic phase. 

We may consider that during the kinetic phase 
the magnetic field intensity increases as B{t) = 
i?oexp (t/rd), being the timescale = ld/6vi^, with 
Svi^ the turbulent amplitude at the viscous scale Id- 

1/2 

The saturation occurs when B{tsat) — Svi^ ( 47 rp) ' ~ 

( 47 rp)^^^, which results in'^: 

Gat ~ In , (I) 

where p stands for the mass density of the fluid, and 
AIa L = (5ul (47rp)^^^/Ho the Alfven Mach number at 
the large scale L with respect to the seed field intensity 
Ho- Large turbulent Reynolds numbers result in faster 
SSDs, however the saturation occurs at lower amplifi¬ 
cation levels reducing its effectiveness. In the opposite 

^ when the magnetic energy density is small compared to the 
kinetic energy density of the smallest turbulent scale 

^ Here we define the Reynolds number as Re = Svi^Lju, being u 
the fluid viscosity and L the largest turbulent scale. We also make 
use of the Kolmogorov scaling relation for the velocity dui oc dU 
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trend, relatively smaller values of Re result in slower 
SSDs, but with larger B{tsat)- However, in order to ex¬ 
plain the cosmological amplification of magnetic fields 
after recombination era we need large B(tsat) and fast 
amplification. During the formation of clusters of galax¬ 
ies, starting at at z ^ 1.0 (which represents a look back 
time of < SGyrs®), velocities of the order of hundreds of 
km/s® are driven at scales as large as L ~ IMpc. For typ¬ 
ical intergalactic gas densities of n ~ 10“^cm“^ one finds 
B{tsat) ~ AtG for i?e ~ 10 — 1000, which resembles the 
values expected for the IGM using the Spitzer viscosity. 
For a seed field Bq .^nG, using the same parameters, one 
also finds Gat > lOGyrs. It is clear, at least from these 
crude estimates, that only superestimated values of the 
magnetic seed field and a fine tunned Reynolds number 
could explain the magneti c field inte nsities observed in 
the local Universe (see e.g. IChdl2014ll . 

Given the inefficient amplification of the field by the 
turbulent dynamo, another mechanism must be con¬ 
sidered. It has been recently pointed that the IGM 
may not be well described by a standard MHD the¬ 
ory. The intergalactic plasma, extremely rarefied, be¬ 
haves as a gyrotropic collisionless plasma, i.e. the Lar- 
mor radius of the ions (ru eB/mcvni^i) is smaller 
than the mean free path Amfp ~ (ngjj)^/^, being g,, the 
ion-ion collision cross-se ction (ISchekochihin et ^ 120051 : 
iSantos-Lima et ^l2014f) . This relation can be rewritten 
as the condition for gyrotropic plasmas as i3gyro(G) > 
10“^°n(cm“^)T(eV)“^/^. In the standard ACDM Uni¬ 
verse, with an average density of baryons at z ~ 1000 
of n^iooo ~ 10®cm“^ and kTziooo ~ leV, one finds ru < 
Amfp for magnetic fields larger than ~ 10“^^G, which 
decreases with time as the Universe expands. Therefore, 
for the given seed fields above, it is natural to assume 
the intergalactic plasma as gyrotropic. If not during its 
whole evolution post-recombination evolution, at least 
from the time when the condition ru < Amfn is satisfied . 

Even in the fluid approximation (iGhew et al.l I1956L 
CGL-MHD hereafter) gyrotropic plasmas are subject 
to a number of instabilities that provide fast in- 
terchange of internal, k inetic and magnetic energies. 
ISchekochihin fc CowlevI (|2006a[l presented a simplified 
analytical study of the time evolution of the magnetic 
energy in the IGM based on such instabilities. Their 
conclusion was that collisionless plasma instabilities may 
result in an “explosive” growth of magnetic energies in 
very short timescales compared to the cosmological evo¬ 
lution of the system. However, such idea has not yet 
been extended into a full theory nor tested numerically, 
which is the main goal of this work. 

In the present paper we study numerically the am¬ 
plification of magnetic fields in turbulent collisionless 
plasmas, under the CGL-MHD formalism. In order to 
achieve this goal we provide a number of direct numeri¬ 
cal simulations of plasma flows, with forced turbulence, 
threaded initially by extremely weak seed fields. Colli¬ 
sionless plasma instabilities are included. The basic the¬ 
oretical aspects of such problem is discussed in Section 2. 

® For a standard ACDM cosmological model, assuming a dark 
energy density A = 0.714, matter density (baryons + dark matter) 
Om = 0.286 and a Hubble parameter h = 0.7 

® as estimated from equipartition with thermal component, i.e. 
mildly transonic turbulence 


The numerical setup and algorithms used are provided 
in Section 3. In Section 4 we present the main results, 
followed by the Discussions and the Conclusions. 

2. MAGNETOHYDRODYNAMIC DESCRIPTION OF 
COLLIONLESS PLASMAS 

The proper description of collisionless plasmas relies 
on the full calculation of the particle dynamics, includ¬ 
ing both eletromagnetic and collisional forces. Unfor¬ 
tunately such approach is of little practical use, with no 
analytical model available yet. One of the alternatives for 
this problem relies on numerically integrating the equa¬ 
tion of motion of an ensemble of charged test particles 
(ions), coupled to a fluid (electrons). Such approach is 
also called particle-in-cell (PIC) numerical simulations. 
PIC simulations provide the dynamical evolution of mo¬ 
mentum distributions of particles as they interact with 
magnetic fields and due to collisions, which allow us to 
directly study e.g. the effects of collisionless plasma in¬ 
stabilities on the isotropization of pressures, or on the 
rise of pressure anisotropies. However, the spatial cov¬ 
erage of the computational domain in PIC simulations 
is limited to a finite number of Larmor radii. For this 
reason it is not possible to study both the evolution of 
the distributions of particle momenta and the system dy¬ 
namics at large scales, such in the case of the IGM, simul¬ 
taneously. At first approximation, we may consider the 
distribution of momenta of particles to be Maxwellian 
(or bi-Maxwellian in the case of magnetized collision¬ 
less plasmas), for which a fluid description of a pr essure- 
anisotropic plasma is available (|Chew et al.lll956D . 

2 .1. Single fluid approximation of a plasma with 
pressure anisotropy 

The derivation of the CGL-MHD equations from the 
Vlasov-Maxwel l equations is provided, for example, in 
iKulsrudI dull). By neglecting heat conduction and 
other possible heating/cooling sources, the one fluid 
CGL-MHD set of equations for a plasma with pressure 
anistropy can be rewritten, in conservative form, as: 


‘ p' 


pn 

pu 

+ V- 

puu -1- Hp -1 - Hb 

B 


Bu- uB 

e 


_eu -1- u • (Hp -1 - Hb). 


■ 0 ■ 
f 


f • V 

where p, u, B, p_L || represent the plasma density, veloc¬ 
ity, magnetic field, and perpendicular/parallel pressures 
with respect to orientation of the local magnetic field, 
respectively, and e = (p_L + p\\/2 -|- pit^/2 -|- B‘^f2) is to¬ 
tal energy density. The anisotropy in pressure for a gy¬ 
rotropic plasma is then defined as A = (j)±—p\\)/p±. Hp 
and Hb are the pressure and the magnetic stress tensors, 
respectively, given as: 

Hp = PJ.I-h (p|| - p_L)bb, (3) 

and 


Hb = (HVStt)! - BB/47 r, 


(4) 
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where I is the unitary dyadic tensor and b = ’RjB. The 
source term f, in the equations above, represents the ex¬ 
ternal force responsible for driving the turbulence. 

The set of equations given in Eql2]is not closed. Since 
the perpendicular and parallel pressures can evolve dif¬ 
ferently from each other, another equation relating the 
time evolution of p_L and p|| is needed. The origi¬ 
nal do uble-adiabatic condition proposed in IChew et al.l 
(1195611 . based on the conservation of first and second 
magnetic invariants, leads to: 


threshold of the actual mirror instability obtained from 
the kinetic theory. This offset is a well-known issue in 
the CGL-MHD closure, but this is not relevant for this 
work as we focus on the evolution of the firehose mode 
only. 

3. AN APPROXIMATE ANALYTICAL SOLUTION FOR THE 
AMPLIFICATION OF B IN THE UNSTABLE REGIME 

By differentiating the pressure anisotropy A in time 
one obtains: 


d 

dt 



= 0 , 


dt\p^) 


= 0 . 


(5) 


1 L dp^ dp\\ 

dt p± dt dt 


(9) 


The CGL-MHD double-adiabatic closure equations 
above may not be exactly correct. It is known that sev¬ 
eral processes, e.g. collisions or enhanced p article scat¬ 
tering by magnetic mirrors (or firehoses le.g. lKunz et al.l 
1201411 1. break these invariants. Another possible closure 
for Eql2]is obtained by simply assuming that the system 
reaches a pressure anisotropy equilibrium fast enough, 
so A = const, during the whole evolution. The closure 
problem is discussed in more details in the following sec¬ 
tions. 


The first term of the right side of Eql5]can be evaluated 
in terms of the magnetic field intensity. The magnetic 
moment, or first adiabatic invariant, p = rui vl/2B = 
constant, of a charged particle subject to the magnetic 
force, is also quasi-invariant for a plasma. Therefore, 
by approximating v to the thermal velocity Uth, associ¬ 
ated to the Maxwellian distribution of velocity compo¬ 
nents perpendicular to the magnetic field, and assuming 
dp/p <C dB/B, we find dp±/p± ~ dB/B, and EqlH] re¬ 
sults in: 


2.2. Wave modes and stability eondition in gyrotropie 

plasmas 

A linear perturbation analysis of Eql2] results in 
the known modified dispersion re lations in gyr o tropie 
plasmas given b elow (see e.g. iHa.n fc Wan^ [2nn^ 
iKowal et al.ll20HT) : 



( 6 ) 


corresponding to the Alfven transversal mode, with 
cosO = k • 'R/{kB) (being k the wavevector of the per¬ 
turbation), and: 


f- 

\k 


2 6 ± 1/52 _ 

f.s 2 


(7) 


corresponding the fast (“-I-”) and slow (“-”) magne- 
tosonic waves, being 


b = 


47rp ' p 
2 


2pj. , /2p|| p_L 


— — ) cos^ 0, 
P P 


c = — cos"^ 9 


3P|| 


3P|I: 


—^ I cos 9 - -b -f ( — I sin' 


P±- 


The equations above result in the following stability 
conditions: 


P± 

for (p||>p±). 


and 

i^>l + /3l\ 

6p|| 

for (p|| < p_l) (8) 


which correspond to the firehose and fluid mirror insta¬ 
bilities, respectively. We must point out that there is a 
factor of 1/6 for the later, which is not present in the 


where 


dA 

dt 


, ^, 1 dB 




1 dpii 
pj_ dt 


Di|A| 


( 10 ) 


( 11 ) 


represents an effective isotropization rate {veS 
J^isoA”^). Notice that the first term in the right hand 
side of the equation above corresponds to the time evo¬ 
lution of p|| as a consequence of the changes in B and 
p^. From here we describe the interchange of parallel 
and perpendicular pressures, mediated by magnetic fluc¬ 
tuations, as the net magnetic scattering of the particle 
distributions (zzscatt P'j}^~^dp\\/dt). The last term 
of the equation above has been deliberately added to ac¬ 
count for the effect of collisions, by means of a Braginskii 
collision frequency va. From these two terms it is possible 
to consider two different regimes with: i) collision dom¬ 
inated isotropization, for ^ r'seatt; and ii) scattering 
dominated isotropization, for <C t'seatt- 
The effect of collisionless instabilities on the pitch angle 
(magnetic) scattering is not fully understood. If we take 
the firehose instability as example, perturbations (e.g. in 
the velocity field) may result in local growth of magnetic 
held intensity (SB). This in turn results in an increase of 
the perpendicular velocities of the ions, i.e. an increase 
of p±, at the expense of a decrease of p||. This could in 
principle be used to estimate dp\\/dt^. However the to¬ 
tal internal energy may not be conserved if particles are 
lost during scattering, or if other energy loss processes 
are included, making the problem not practical without 
a f ull kinetic description. An alt ernative, as pointed 
by ISchekochihin fc CowlevI (l2006bll . is to assume this 
term to saturate the ampliheation of the magnetic held 
perturbations quasi-linearly at the wavelength of fastest 
growth, which occurs around the Larmor radius (w ^ Oi, 


^ e.g. as with the use of the 2nd adiabatic invariant 
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being flrmi = eB/m-iC the cyclotron frequency of the 
ions ) (|Garv et al.lll998ll . therefore 5B/B\ sat ; 

given the maximum growth rate (see Eqsiniand[7]): 

7max (^|A| - Gi, (12) 

with Of = 1 and 1/2 for mirror and firehose instabilities, 
respectively. As we show further in this paper the fire¬ 
hose instability is of major interest in the amplification 
of B^ and we use a = 1/2 from here. The scattering fre¬ 
quency can then be estimated as I'scatt — SB'^ / 

(|A| — 2//3)^/^fli, and Eg ITOl becomes: 

3.1. Magnetic field amplification with constant A 

This problem is tractable analytically if we assume 
that the evolution of the firehose unstable regions in the 
plasma occur with constant A. This approximation is 
justified by the fact that any process that may amplify 
the pressure anisotropy (e.g. turbulence, anisotropic cos¬ 
mic ray scattering, and others) is precisely counterbal¬ 
anced by the ones leading to isotropization. This is ac¬ 
tually a probable scenario given that the increase of the 
pressure anisotropy by external sources is independent 
on A, while the isotropization processes are a function of 
the anisotropy itself. Therefore, the system should evolve 
around an “equilibrium” anisotropy value Aq. If the re¬ 
laxation timescale is short compared to the dynamical 
timescales of the system, we may assume A(t) ~ Aq. We 
must stress that the validity of such conjecture depends 
on the level of anisotropy itself, i.e. it is probably correct 
for Aq —>■ 0 and certainly does not stand for Aq ^ 0. 


smallest turbulent scales. However, as B grows, —>■ 

and the plasma becomes gyrotropic. This is supposed 
to occur for extremely weak magnetic fields (i?gyro(G) > 
10“^°n(cm“^)r(eV)“^/^), possibly even earlier than the 
saturation of the turbulent dynamo at small scales. 

Another interesting aspect of the turbulence during 
the early evolutionary phase is to generate the pressure 
anisotropy by itself (EafT^. This is physically under¬ 
stood as the consequence of a fast growth of the mag¬ 
netic field not counterbalanced by collisions. Depending 
on the properties of the turbulence, one obtains | A| ^ 0. 
Such values are not expected in highly magnetized plas¬ 
mas. Once Vii r'scatt, due to the amplification of 
B, the balance is reached since r'scatt strongly depends 
on the magnetic field intensity. Local measurements of 
the Earth’s magnetosphere and the solar wind reveal 
pressure anisotropies around its quasi-stability threshold 
|A| ^ |A|crit, and many signatures of its role on the dy- 
namical evolution of these systems have been found (e.g. 
Winterhalter et all 119951 : ISoucek et al.ll200§ : iBale et 
2009t) . The interesting implication of Eqlin] is that, due 
to its lower magnetization level (much higher /3), |A| in 
the IGM/IGM could have been much larger at higher 
redshits, compared to the values found presently in the 
Solar System. 


3.1.2. the case of va < r^scatt 

Once the magnetic field becomes strong enough, for 
the plasma to become gyrotropic, instabilities will take 
over the amplification of the magnetic field. Now, as 
A = const and <C Vscatt, and replacing the cyclotron 
frequency Eg fOl becomes^: 


dB e 

dt TOiC 1 -I- |A| 


(16) 


3.1.1. the case of va > Vacatt 

Let us first consider the case in which collisions domi¬ 
nate the isotropization of pressure. This limit could well 
describe the initial stages of the evolution of a cosmo¬ 
logical seed field after the recombination era, given that 
(3 :$> 1 and va Gj. The evolution of the magnetic 
field in this case depends mostly on the dynamics of the 
plasma, and not on instabilities driven by the pressure 
anisotropy. Given that the resistive dissipation is ex¬ 
tremely small in the intergalactic medium we may ap¬ 
ply the “frozen-in” condition to the plasma. The evo¬ 
lution of the magnetic field is then dominated by the 
turbulent rate of strain at th e viscous scales (see e.g. 
ISchekochihin fc Gowlevll^OOfibll : 

— ^ ~bb-Vu-V-u- —(14) 
B dt L 

For A = const and va ^ r'scatt, and combined with 
EofTH EqllSlis then reduced to: 



As discussed in the previous section the kinetic phase 
of the turbulent dynamo fEo lldl) provides an exponen¬ 
tial growth of the seed field up to equipartition at the 


which, at fi ^ |A|~^ limit, result in a first order non¬ 
linear ordinary differential equation with an analytical 
solution: 

Bit) {B^^ - At)-\ (17) 

I 13/2 

with A ~ i-HA| • magnetic field growth is ex¬ 

plosive around t ~ A~^, and the saturation occurs for 
B(rgr) ^ Bq. The quasi-stability condition, given by 
(3 ~ 2/1 A| level, is reached in a timescale: 

rg,(s) ^ (|A|’>G,.o)-i - 10-X'(G)|A|-’^, (18) 

where represents the cyclotron frequency measured 
for the seed field Hq, and = 1/2, or 3/2, for |A| 1 

and |A| ^ 1, respectively. EolTSl then provides an ap¬ 
proximate timescale for the explosive growth of the mag¬ 
netic field due to the collisionless plasma instability. For 
example, if we consider a seed field ofi?o ~ 10“^^G sub¬ 
ject to such instabilities, with |A| > 10“^, one finds am¬ 
plification up to equipartition with pressure anisotropy 
(/3 ^ |2A|~^) at t < 400Myrs. 

® notice that the amplification of B in our model can occur for 
firehose instability, i.e. py > (see Eql9]l, and therefore A < 0 
and 1 — A > 1 
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Fig. 1. — Time evo luti on of magnetic field a mplitu de in the 
approximate case of Ea ll6l as described in Section l3.1.2l 

As pointed before, such values of pressure anisotropies 
are easily reached in turbulent weakly magnetized col¬ 
lisionless plasmas. The main conclusions of this sec¬ 
tion are that: i) pressure anisotropies could be naturally 
generated in the IGM/ICM in the early Universe (post¬ 
recombination era) at timescales comparable to the eddy 
turnover time at viscous dissipation scales, and ii) the in¬ 
stabilities driven by the pressure anisotropy . A full nu¬ 
merical integration of Eallbl is shown in Figure [TJ which 
confirms the main results obtained in the approximate 
limits /3 » |A|~^ and S(Tgr) ^ Bq. 

4. NUMERICAL SIMULATIONS 

In the previous section we demonstrated that the 
evolution of plasma instabilities driven by a pressure 
anisotropy could, in principle, explain the amplification 
of a magnetic field seed at recombination era to the val¬ 
ues observed in the local Universe. These estimates, 
however, were obtained under strict approximations and 
limit cases. In this section we describe the numerical 
methods used for modelling of the dynamical evolution 
of the unstable plasmas. 

4.1. Governing equations and methods 

In order to perform the numerical modeling of the 
plasma evolution with the constant pressure anisotropy 
and turbulence injection we used the GODUNOV code® 
which solves the CGL-MHD equations (Eqs. [2|) in the 
conserved form. 

For simplicity py = and p± = a\p are expressed by 
the sound speeds ay and ox, parallel and perpendicular 
to B, respectively. Therefore the momentum equation 
can be rewritten as 


dpvL 

TT 


-b V' 


puu -b 


a\p' 




= f, 


(19) 

where a = (a| - a\)/Vl = \ (/3|| - /3x) = i/3x (^ - 1) is 
the pressure anisotropy degree with ^ = py/px = «y/ax 
being the pressure ratio. 


^ The MHD version of the code is publicly available from the 
website http: //amuncode. org 


The numerical integration of the modified^® CGL- 
MHD equations, was done in GODUNOV using 
the second-order s hock-capturing Godunov-scheme (see 
iKowal et al.ll20Illl . The time integration was done us¬ 
ing the second order Strong St ability Preserving Ru nge- 
Kutta (SSPRK) method (see iGottiieb et al.l l2009l and 
references therein). Spatial reconstruction was done 
using the second order the total variation di minishing 
(TVD ) interpolation with the Van Leer limiter (|van Leer! 
Il974f) . and numerical fluxes were calculated using the 
general Harten -Lax-van Leer (HLL) Riemann solver (see 
lEinfeldtlll988l e.g.). We incorporated the fi eld interpo- 
lated constrained transport (CT) scheme lsee lT6thll200fll 
e.g.) into the integration of the induction equation to 
maintain the V • B = 0 constraint numerically. 

4.2. Model of turbulence 

In our numerical mo deling we drive turbulence using 
a method described bv lAlveliui (II999D . The forcing is 
implemented in spectral space and concentrated around 
the injection scale related to a wave vector kinj. We 
perturb a finite number of discrete Fourier components 
of velocity in a shell extending from kinj — ^kinj to kinj + 
Akinj with a Gaussian profile of the half width kc and the 
peak amplitude u/ at the injection scale. The amplitude 
of driving is solely determined by its power Pinj- The 
parameters describing our forcing do not change during 
the evolution of the system. 

The randomness in time makes the force neutral in the 
sense that it does not directly correlate with any of the 
time scales of the turbulent flow, and it also determines 
the power input solely by the force-force correlation. In 
the models presented in this paper we use isotropic forc- 
ing. 

In particular, the total amount of power input from 
the forcing can be set to balance a desired dissipation at 
a statistically stationary state. In order to contribute to 
the input power in the discrete equations from the force- 
force correlation only, the force is determined so that the 
velocity-force correlation vanishes for each Fourier mode. 
The procedure of reduc i ng th e velocity-force correlation 
is described in lAlveliusI (|1999ll . 

In Eq. (froi) . the forcing is represented by a function f = 
pa, where p is local density and a is random acceleration 
calculated using the method described above. 

4.3. Initial and boundary conditions 

For our calculat ions, similar to our earlier studies (see 
iKowal et al.ll20li]l . the initial pressures and the uniform 
initial field Bq are the controlling parameters. We define 
the Alfvenic Mach number of the injected turbulence as 
Ma = {5v/ca)- The angle brackets () represent the 
volume average of the parameter. 

The distribution of p = 1.0, v = 0 and B = Bqx, is 
uniform in the whole domain. We do not set the viscosity 
nor the resistivity coefficients explicitly in our models, 
and dissipation is determined by the numerical scheme 
only. Therefore the scales at which the dissipation starts 
to be important is defined by the numerical diffusivity of 
the scheme. 

standard solution of CGL-MHD equations envolves the use of 
the first and second adiabatic invariants as closures for the prob¬ 
lem. Here, this is slightly modified as we avoid the use of the 2nd 
adiabatic invariant by fixing the pressure anisotropy as constant. 
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TABLE 1 

Initial and saturation properties of the numerical simulations 


model 

res 

Ma,o “ 

00 

P± 

P|| 

|A|‘ 


"Tsat (^dyn) 

fh512Ma0D1.0 

512^ 

2.0 

2.0 

0.25 

1.0 

3.0 

1.6 

0.52 

mr512Ma0D0.5 

5123 

1.0 

3.0 

4.0 

1.0 

0.75 

4.3 

- 

mr256Ma4D0.8 

2503 

10"* 

6.8 X lO'^ 

1.0 

0.04 

0.96 

2.5 X lO'^ 

- 

mr256Ma4D0.6 

2503 

lO'’ 

7.5 X 10'^ 

1.0 

0.25 

0.75 

2.2 X 10'^ 

- 

fli512Ma4D1.0 

5123 

lO'^ 

2.0 X 10® 

1.0 

4.0 

3.0 

7.6 

0.36 

fh256Ma4D1.0 

2503 

10^ 

2.0 X 10® 

1.0 

4.0 

3.0 

7.8 

0.62 

fhl28Ma4D1.0 

1283 

lO'^ 

2.0 X 10® 

1.0 

4.0 

3.0 

7.9 

1.40 

fh256Ma4D4.0 

2503 

lO'^ 

9.0 X 10® 

1.0 

25.0 

24.0 

1.4 

0.13 

fh256Ma4D9.0 

2503 

lO'^ 

3.4 X 10® 

1.0 

100.0 

99.0 

0.6 

0.07 

fh256Ma6D4.0 

2503 

10® 

9.0 X IQi® 

1.0 

25.0 

24.0 

1.3 

0.18 

fh256Ma8D4.0 

2503 

10® 

9.0 X IQi® 

1.0 

25.0 

24.0 

1.4 

0.22 


“MIa.o = is the Alfven Mach number of the turbu¬ 

lence injected with respect to the initial magnetic field, i.e at t = 0 


"A = (px -P||)/pt 

The simulation box is perfectly periodic in all direc¬ 
tions. 

4.4. Simulated models 

For the present work we performed a number of numer¬ 
ical simulation with different plasma initial parameters, 
as well as different stability regimes, as given in Table 
[TJ We also performed simulations varying the numerical 
resolutions for the purpose of convergence validation. 

5. RESULTS 

A comparison between the spatial distributions of 
the density and velocity fields in CGL-MHD and stan- 
dard MHD t urbul ent models have been p rovid ed in 
iKowal et al.l (|2011ll and ISantos-Lima et ahl (|2014ll . In 
those works it has been shown that CGL-MHD models, 
depending on the filling factor of unstable regimes and on 
the degree of pressure anisotropy itself, may lead to evi¬ 
dent differences in the statistics of turbulent fields when 
compared to standard MHD turbulence. Those studies 
were performed for relatively strong magnetic fields, i.e. 
for px|A| > (H^), and therefore the statistics of B were 
w eakly dependent on CGL -MHD instabilities (see Fig.6 
in iSantos-Lima et al.ll201^ . Differently to the previous 
works, we focused our modelling on the weakly magne¬ 
tized turbulent regimes (p_l|A| 

5.1. Magnetic field struture and statistics 

As the simulations are initiated, the driving of tur¬ 
bulence results in fluctuations of density, and of the 
velocity and magnetic fields, at different scales. It is 
well known for many decades that firehose and mirror 
instabilities, given the growth rate dependence with k 
(Eqsin]and[7]), result in the emergence of structures much 
smaller than those observed in standard MHD turbulence 
(|Kowal et al.l 1201 ill . In standard MHD, the forcing, i.e. 
the process of injection of turbulence, is the dominant 
dynamical process and therefore dominate the process of 
structure formation. In the collisionless plasma approxi¬ 
mation, on the other hand, the growth rate of the insta¬ 
bilities rises indefinitely with the wavenumber of the fluc¬ 
tuation, up to fcmax ^ Cl. The numerical scheme imple¬ 
mented in this work for solving the CGL-MHD equations 
do not account for finite Larmor radius effects. There¬ 
fore, the maximum growth rate of the instabilities must 
be related to the minimum lenghtscales of the system. 


which is the size of the finite volume of the space dis¬ 
cretization, i.e. the grid cell. This means that both the 
growth rate and the geometry of the amplified magnetic 
fields are resolution dependent. 

In Figure [5] we present the aparent configuration of the 
magnetic field lines in the mid-slice (in the z-direction). 
The geometry of the field lines is illustrated by the line 
integral convolution (LIC) technique. The LIC imaging 
algorithm consists of rendering a map of random stream¬ 
lines that follow the orientation of the local field. The 
values are normalized to arbitrary units in order to pro¬ 
vide a texture map. Such texture map provides a visual 
qualitative measure for the scaling of the magnetic field 
lines. Figure [5] shows the LIC technique applied for the 
the magnetic field lines at saturated stage for |A| = 3.0 
and /3o = 2x 10®, corresponding to firehose unstable mod¬ 
els, for the given resolutions of 128® (left), 256® (center) 
and 512® (right). 

It is visually clear from these plots that the finer the 
resolution the more fiocculent the geometry of the am¬ 
plified magnetic field becomes. This result is expected 
since for |A| = 3.0 and /3o = 2 x 10® models the thermal 
pressure, including the pressure anisotropy free-energy, 
is much larger than the magnetic and turbulent driving 
energy densities. Therefore the firehose instability domi¬ 
nates the plasma dynamics and the structure formation, 
which growth rate should peak around the smallest scales 
available in the system. This is also quantitatively mea¬ 
sured by means of the normalized power spectrum of the 
magnetic field, shown in Figure |3] (top). The three di¬ 
mensional power spectra of the magnetic field for these 
models reveal that the peak shifts towards larger fc’s as 
the numerical resolution increases. The peak shift is lin¬ 
ear with the resolution, ocurring at rounded values of 
kpeak —16, 33 and 65, for resolutions of 128® (black), 
256® (red) and 512® (blue), respectively. Perturbations 
smaller than 6 — 8 grid cells suffer mostly the effects of 
numerical dissipation. 

Notice that, theoretically, the maximum growth rate 
should occur at fc ^ fcpeak but numerical diffusion is re¬ 
sponsible for the dissipation of coherent structures at 
the smallest scales. Both the numerical viscosity and 
resistivity act together to reduce the effect of the in¬ 
stabilities in the growth of velocity and magnetic field 
perturbations. If the numerical viscosity is considered 
to be Laplacian we can, in principle, estimate the peak 
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Fig. 2.— Central slice line integral convolution (LIC) map of the magnetic field lines at saturated stage for |A| = 3.0 and /3o = 2 x 10®, 
firehose unstable, models with different resolutions: 128® (left), 256® (center) and 512® (right). 


location as the lengthscale at which the numerical vis- 
cous loss rate (47r^fcpg^j^z/num) equals the instability grow 

rate, which results in fcpeak y^|A| — 2/3“i(87r^z/num)~^- 
Given that the asymptotic values of /3 depends only on 
A, when saturation of instabilities occurs the peak lo¬ 
cation becomes insensitive to (|A| — 2/3“^), resulting in 
fcpeak oc i-®- approximately linearly proportional to 

the numerical resolution as suggested before. 

Also, it is noticeable that the numerical resolution, ex¬ 
cept for the location of the peak, is not relevant for the 
general shape of the power spectrum of the magnetic field 
fluctuations. A power spectrum P oc with a ~ 2, is 
observed in the range of scales ~ [fcinj, fcpeak], in agree¬ 
ment with the saturation estimates for 6B^ oc 7 I ex 
(see Section 3), for all the numerical resolutions. This is 
a good indicative of convergence for the models. The ob¬ 
tained power spectrum is steeper than the one expected 
for the kinetic-phase of the small-scale dynamo, which 
shows that the instabilities are the dominant process in 
the field amplification here. This also show that the mag¬ 
netic field at large scale is comparatively weaker. 

The effects of | A| on the structure of the magnetic field 
is also shown in Figure [3] (center). Here we present the 
spatial spectral power of the magnetic field for the 5 dif¬ 
ferent models with same thermal perpendicular pressure 
(p_l) and numerical resolution (256^), at t = 2.0tdyn, but 
for different parallel pressures (py), resulting in mirror 
(“M”) and firehose (“FH”) unstable systems with differ¬ 
ent anisotropies. The mirror unstable models, in contrast 
to the firehose unstable ones, present a decreasing power 
with wavenumber. This is because, as shown above, the 
amplification of the magnetic fluctuations in the firehose 
regime is dominated by the instabilities (positive slope 
of the power spectra), while in the mirror unstable re¬ 
gions it is basically subject to the dynamics of the flow, 
i.e. the turbulence (responsible for the negative slope 
of the power spectra). Also, as we discuss in more de¬ 
tails in the next subsection, the amplitude of the power 
spectrum is a function of the pressure anisotropy for the 
firehose unstable models, while it is insensitive to |A| 
in the mirror unstable models. In the firehose case, the 
larger the anisotropy the larger is the total power of the 
magnetic fluctuations, which is related to the real val¬ 
ues of the fluctuations and to the degree of ampliheation 
itself. 


Since saturation is an important condition for the iden¬ 
tification of the peak and total power of the magnetic 
fluctuations we studied the time evolution of the spa¬ 
tial power spectrum of B of one of the firehose unstable 
models (|A| = 24,/3o = 9 x 10^® and 256^). We chose 
the model with the smallest magnetic seed field given its 
longer saturation timescale, which could eventually in¬ 
fluence the spectral power distribution. As can be seen 
from Figure [3] (bottom), the spectra are very similar for 
all snapshots {t = 0.5,1.0,1.5 and 2.0tdyn), with exactly 
the same slopes and peak locations. This means that 
the saturation occurs at t <g; fdyn, and that the turbu¬ 
lence driven - even though super-Alfvenic (with respect 
to the seed field) - is not able to modify the spatial dis¬ 
tribution of the amplified field fluctuations. It is worth 
mentioning here that, since the non-linear regime (as for 
the small-scale dynamo) apparently does not exist, an¬ 
other mechanism is therefore required to continue the 
magnetic field amplification on larger scales. 

5.2. Magnetic field amplification 

As indicated by the power spectra shown in Figl2 
the firehose unstable models result in amplihed magnetic 
field fluctuations. In the left panel of Fig0] we present 
time evolution of the magnetic energy density for mod¬ 
els with different pressure anisotropy ratios. The mag¬ 
netic energy density was averaged over all grid cells of the 
simulated box. As mentioned above, the mirror unsta¬ 
ble models show moderate growth rates for the magnetic 
field perturbations, which occur at timescales of ~ tdyn, 
in agreement with a typical slow turbulent dynamo. The 
firehose unstable models on the other hand present much 
larger amplification values, corresponding to > 7 orders 
of magnitude for the three models shown. Also, the am¬ 
plification timescales are much shorter than in the mirror 
unstable case, with Tgr <C tdyn- For these models the sat¬ 
uration of magnetic energy density occurs at larger levels 
for larger pressure anisotropies, as expected for the satu¬ 
ration threshold /3“J oc A| (see Ea fT6)) . Also, in qualita¬ 
tive agreement with Eo llSl the saturation timescales for 
the magnetic field amplification decreases with pressure 
anisotropy. The timescales and amplification levels for 
saturation are also presented in Table [TJ 

Another interesting feature related to the time evolu¬ 
tion of the magnetic energy density is the presence of 
a “knee” separating two different growth regimes. The 
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Fig. 3.— Spatial power spectra of the magnetic field for: different 
numerical resolutions (top), different pressure anisotropies (center), 
and at different stages of the system evolution (bottom). 

energy density level of the knee is similar for all mod¬ 
els. As pointed in Section 3, if the turbulent rate of 
strain is larger than the growth rate of the instabilities 
the turbulent dynamo regime dominates. This is pre¬ 
cisely the physical motivation of the knee. At the early 
stages turbulence is driven at large scales, which are re¬ 
lated to small growth rates of the instabilities. At this 
point Eq[T4] is not yet valid given that the turbulence is 
not fully evolved, but could be replaced with a reduced 
effective Re number instead due to the shorter inertial 
range. The actual magnetic growth rate, related to the 
rate of strain, increases with time as the cascade devel¬ 
ops and the inertial range grows. This behavior is also 


observed in the curves of magnetic energy density evolu¬ 
tion. Once perturbations grow at a particular scale, at 
which the growth rate of the instabilities take over the 
numerical dissipation, the explosive amplification takes 
place. 

The right panel of Fig|T] presents the time evolution 
of the turbulent Alfvenic Mach number for the differ¬ 
ent models. All models start with similar turbulent 
strengths. From this it is clear that, for the cases with 
A < 0, the bulk of the magnetic field amplification is not 
related to the turbulence itself (though it is needed to 
seed the fluctuations that are subject to the instabilities 
at a later stage). The mirror unstable models demon¬ 
strate very similar evolutions while, on the other hand, 
the firehose unstable models show different evolution de¬ 
pending on the level of pressure anisotropy. Again, as 
pointed before, the larger the pressure anisotropy the 
faster the plasma reaches saturation. 

The same study is presented in Fig. [51 but for the dif¬ 
ferent numerical resolutions. We can observe that the 
evolution of the average magnetic energy density and 
the Alfvenic Mach number present similar profiles for the 
different models. Therefore, the overall time dependency 
function of these quantities is well converged in our mod¬ 
els. The convergence is good for the asymptotic value 
of the average magnetic energy density as well, which 
saturates for similar values regarless the resolution and 
depends on A exclusively. The timescales for satura¬ 
tion however is resolution dependent, as we have pointed 
previously, being shorter for higher resolutions. This is 
in agreement to the fact that finer resolutions allow the 
growth of the shorter wavelength perturbations, which 
have faster growth rates. 

The effect of the seed field intensity on the mag¬ 
netic growth rate is shown in FigjSl for the models with 
A = 24.0 and numerical resolution of 256^ cells. Here the 
different seed field intensities correspond to initial ther¬ 
mal to magnetic pressure ratios f3o = 9 x 10^^ (red) and 
/3o = 9 X 10® (black). The time evolution of the Alfvenic 
Mach number (right panel) shows that saturation occurs 
at similar levels regardless of the seed field intensity. The 
saturation occurs once instabilities drive the amplifica¬ 
tion to ^ |A|pj^ level, as predicted from the CGL-MHD 
stability condition. Weaker initial magnetic fields result 
in later saturation timescales. The delay, compared to 
the models with stronger initial fields, occurs as the road 
up to the saturation level must take longer for weaker 
seed fields. 

Notice that the “knees”, i.e. the transition between 
the turbulent and instability driven dynamos, occur at 
similar times (t ~ 0.02 — O.OStdyn) regardless the initial 
magnetic field intensity. Also, the magnetic field ampli¬ 
fication, relative to the seed field intensity, of the knee is 
similar in the three models. Both features are in agree¬ 
ment to FqHTl That equation reveals that, for the same 
timescales and similar turbulent properties, the relative 
amplification AB/B should be similar for all models. 

The time evolution of the minimum and maximum 
magnetic energy densities are also shown in the left panel 
of FigjOl as dashed lines. These follow qualitatively the 
same profiles of the mean magnetic energy density, ex¬ 
cept for the short period before the “knee”, when the 
maximum energy increases while the minimum energy 
decreases with time. This process is due to the tur- 
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Fig. 4. — The time evolution of the magnetic energy density normalized by its initial value (left) and of the Alfvenic Mach number 
(right). The data shown correspond to the models with px = 1.0, 256® resolution, but py = 100.0 (blue), 25.0 (green), 4.0 (yellow), 0.25 
(red) and 0.04 (black). 




Fig. 5. — Same as Figure|4]but for the models with A = 3.0, /3o = 2 X 10®, and different resolutions: 128® (blue), 256® (red) and 512® 
(black). 




Fig. 6.— Same as Figure|4]but for the models with A = 24.0, 256® resolution, and different seed fields: /3o = 7 X 10^® (blue), /3o = 9 X 10®® 
(red) and /Ig = 9 X 10® (black). 
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bulence that increases the dispersion of the magnetic 
field distribution, before the growth of the instabilities. 
The second phase happens when the instability growth 
rate becomes more important than the turbulent rate of 
strain, therefore the “knee” occurs earlier for the maxi¬ 
mum value curve (larger rate of strain), and later for the 
minimum value curve (lower rate of strain). 

The turbulent broadening of the magnetic energy den¬ 
sity is clearly seen in the probability distribution func¬ 
tion (PDF) shown in FiglT] In the left panel we present 
the PDFs as they depend on the pressure anisotropy. 
We show the models with pj_ = 1.0, 256^ resolution, for 
P|l = 100.0 (black), 25.0 (red), 4.0 (yellow), 0.25 (green) 
and 0.04 (blue), all at 2.0tdyn. All models start with a 
delta function PDF, at B / Bq = 1.0. The broadening 
observed in the mirror unstable models exemplifies the 
role of the injected turbulence on the distribution of the 
magnetic energy over the simulated domain. There is 
clearly no net amplification of the averaged held inten¬ 
sity though. The hrehose unstable models present PDFs 
that are shifted and skewed. Not only the averaged held 
intensity is amplihed, as discussed above, but these PDFs 
- as a whole - are shifted towards larger intensities, being 
the shift proportional to the pressure anisotropy. 

The skewness of the PDFs arises as their peaks are fur¬ 
ther shifted towards larger intensities, compared to the 
bulk of the distribution. This occurs if the amplihcation 
of the magnetic held is larger/faster for stronger helds. 
Such behavior could be understood as a transient pro¬ 
cess due to the turbulent broadening, in which turbulence 
would naturally build-up strong held in more regions as 
time evolves, but this is not the case. The mirror un¬ 
stable models do not present similar skewness. Also, the 
time evolution of the PDFs do not support such possi¬ 
bility. In the right panel of FiglT] we show the time evo¬ 
lution of the PDF of magnetic held intensity for a single 
model. We have selected the model with 128^ resolution 
because it presents the slower evolution of the distribu¬ 
tion, and the PDFs obtained at different snapshots can 
be compared. Here we show four PDFs calculated at 
t = 0.5, 1.0, 1.5 and 2.0fdyn- The general PDF prohles 
change slightly with time, but is similarly skewed from 
t = 0.5tdyn up to saturation time. The aparent constant 
skewness of these curves shows that the changes in the 
distribution of magnetic held intensity occurs too early 
to be accounted for the turbulence, and should be related 
to the instability itself. This is in agreement with EalTBl 
which predicts faster amplihcation for stronger helds. 

5.3. Comparison to the analytical approximations 

Obviously, any comparison between the simulations 
presented here and the magnetization of the intergalac- 
tic medium needs caution. The numerical simulations 
presented here correspond to full three-dimensional and 
stochastic numerical solutions of the problem. It repre¬ 
sents a good benchmark for the analytical estimates (and 
approximations related) of the simplest zero-dimension, 
constant-A, /3 —>■ oo limit, derived in Section 3. 

We start analysing the dependency of the saturation 
timescale observed in the simulations with the pressure 
anisotropy ratio (A). In the top panel of FigI5|we present 
the values obtained from the simulations, as described 
in Table [1] together with a number of analytical solu¬ 
tions for comparison, as given by Ea ll6l The trian¬ 


gles represent the three models with similar parame¬ 
ters except for the pressure anisotropy, therefore we fit 
the analytical solutions to these data. The squares cor¬ 
respond to the models with different initial conditions 
and are included for the sake of comparison only. Here 
we also consider the fact that the simulated instabili¬ 
ties are slightly delayed by the fact that we do not ini¬ 
tiate the models with fully developed turbulence. On 
the contrary, the initial configuration is that of an uni¬ 
form magnetic field, which is disturbed as turbulence 
injected at large scale cascades. This turbulent delay 
(^tturb) is introduced as a free constant parameter, i.e. 
Tgr ~ tsat — ^tturb- Overall there is reasonable agreement 
between the simulations and the analytical solutions, be¬ 
ing those with O.Oltdyn < i^fturb < 0.05tdyn particularly 
good. It is interesting to note that this range is similar 
to the timescales of the “knees” observed in FigUj 
In the lower panel of Figl8|the same study is presented 
but now for the initial values of the seed fields. Here 
the triangles represent the three models with similar ini¬ 
tial conditions except for the seed field (/3o = 9 x 10®, 
9 X 10^^, and 9 x 10^®). According to Eq[T8l the best fit 
to the data should occur for Tgr oc Bq, with 'd = —1.0. 
However, as clearly seen from the plot, the best fits occur 
for a much weaker dependence (or even no-dependence) 
of Tgr with the initial seed field, being -d ~ —0.05. The 
apparent discrepancy between the two results is not real 
though. Notice that EqlTHI is derived as the magnetic 
field evolution once the condition va <C r'scatt is fulfilled. 
The term Bq is therefore somewhat misleading in the 
sense that one must consider the magnetic field intensity 
once z/ii <g; zzscatt, i.e. a transition field intensity. For the 
analytical model the transition field intensity is deter¬ 
mined by z/ii, while in the simulations by the numerical 
dissipation, i.e. in both cases it is insensitive to initial 
magnetic field. Another major effect resposible for this 
discrepancy is that fact that the growthrate depends on 
the Larmor radius, which cannot be propely modelled in 
the simulations. Naturally, Tgr obtained from the simu¬ 
lations is overestimated as B increases with time. 
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Fig. 7.— Probability distribution functions (PDFs) of the magnetic field intensity, normalized by its initial value, computed from all grid 
cells. Left: for the models with different pressure anisotropies, at 2.0£dyn- The data represent the models with p± = 1.0, 256^ resolution, 
but p|| = 100.0 (black), 25.0 (red), 4.0 (yellow), 0.25 (green) and 0.04 (blue). Right: for the = 2 x 10®, A = 3.0 and 128® resolution 
model at the different times: t = 0.5, 1.0, 1.5 and 2.0£dyn- 


6. DISCUSSION 

Several authors have alread y discussed the ef- 
fects of turbulence (see e.g. iBa neriee fc .Tedamzi^ 
2004 ISaveliey et lill 120121: lapichino fc BriiggenI 20121 : 
Saveliev et al.l 120131 : IChol 2014 and many others) 
and structure fo rmation (see e.g. iDolag et al.l 120021 : 
iVazza et al.l 12014 and many others) on the evolution 
of the intergalactic magnetic field. The lack of detec¬ 
tion of TeV radiation from blazars points towards the 
existence of amplified magnetic fields, with filling fac- 
tors around 60%, even at the intergalactic medium voids 
(|Neronov fc Vovldl2010t [Polag et ai.ll2011[l . For this rea¬ 
son, turbulence - as a ubiquitous phenomenon - has 
been preferred as main mechanism for field amplifica¬ 
tion. The common problem related to turbulent models 
of collisonal plasmas is the typical timescale needed for 
magnetic field amplification. As explained previously the 
timescales associated to turbulent amplifications are too 
large for the magnetic seed amplitudes suggested so far. 
Structure formation, on the other hand, would be domi¬ 
nant (due to shear, compression or other star-formation 
related mechanisms) in further amplification of the fields. 

In this work we explore a possible scenario for 
magnetic field amplification at high redshits (post¬ 
recombination) , based on the growth of collision¬ 
less plasma instabilities. Two similar works have 
dealt with the proble m in a similar approach, be- 
ing one analytical (ISchekochihin fc Cowled l2006all 
(see also ISchekochihin fc CowlevI 2006bll . and the 
other purely numerical (ISantos-Lima et ^I20l4) . In 
(ISchekochihin fc CowlevI l2006a^ an analytical simplifi¬ 
cation of the problem is posed, which inspired most of 
the analytical derivation presented here as well. In their 
work, however, it is assumed that both, the damping at 
small scales and the pressure anisotropy A, are similarly 
related to t'seat and i/a. This in order to obtain the time 
evolution of A and B in terms of an effective collision 
rate. In the present we do not constrain z/gcat which is, 
in another hand, implicitly constrained to the assump¬ 
tion that A is kept constant during the evolution of the 
system. Interestingly, for the super-exponential phase, 
both works predict the magnetic field intensity growth 
rate as oc (I — tjtc)~^, though with different slopes and 


characteristic timescales. In both estimates the growth 
of the magnetic field by collisionless instabilities could, in 
principle, explain the magnetization of the local Universe 
based on very weak seeds (< 10“^’^G). Only the analyt¬ 
ical model presented in this work has been compared to 
numerical simulations. 

With res p ect t o the numerical solutions, in 
iKowal et al.l ()2011ll the authors performed simula¬ 
tions with strongly magnetized plasmas and could not 
address the ampliheation problem. In this sense, the 
numerical study presented here stands as an extension of 
that particular work but focused on the statistical prop¬ 
erties of the magnetic field and the role of collisionless 
plasma instabilities on the amplification of the magnetic 
field. The scalings found for density and velocity fields 
in the firehose unstable models of that work resemble 
those obtained here for the magnetic held, with a power 
spectrum peak at the smallest scales available in the 
system. Similar st a tistica l properties were obtained by 
ISantos-Lima et ahl ()20l4) . who performed collisionless 
plasma simulation with variable pressure anisotropy 
and studied the statistical properties of the magnetic 
helds as well. The plasma properties were studied in 
a variable-A framework though, with isotropization 
mechanisms being introduced ad hoc. Still, it was shown 
that even in the cases where isotropization is strong 
enough to wash out the imprints of instabilities in the 
statistical properties of the plasma, their dynamical 
effects are important locally. The ampliheation of the 
magnetic held was studied and, as a consequence of 
small volume coverage of unstable regions, their general 
conclusion was that the instabilities were not efficient. 
Differently to the assumption made in this work, their 
models do not include an external source term to excite 
the pressure anisotropy, which could - in principle - 
change the system to a quasi constant-A case, similar to 
ours. 

Possible observational evidence (or constraints) for the 
process described in this work could be obtained from the 
cosmic microwave background. If collisionality condition 
was satished at pos-equipartition era, the magnetic held 
at small scales could have been amplihed prior to the 
recombination era, and the ionization structure of the 






























































Magnetic field amplification in collisionless plasmas 


13 




Fig. 8. — Comparison between the numerical simulations and the 
analytical predictions given in Section 3. Top: Tgr vs A correlation. 
Triangles correspond to the simulations with similar initial condi¬ 
tions except for the pressure anisotropy (A = 3, 24 and 99), which 
are the ones used for t he correlation analysis. The lines correspond 
to the solution of Ea ll6l for different turbulent delays <5tturb (see 
text). Bottom: Tgp vs Bq correlation. Again, triangle represent the 
models with similar initial conditions, but now with different seed 
field intensities (/3o = 9 X 10®, 9 X 10^^, and 9 X 10^®). The lines 
correspond to correlations Tgr oc Bq , being Eo llSI the case with 

= - 1 . 0 . 

Uni verse would be different, with effects on Silk damping 
(see I Jedamzik &: Abell [2M^ . On the other hand, if the 
plasma becomes gyrotropic at z < Zreci such effect would 
not be relevant. 

The main problem with the models mentioned above is 
the lack of a proper physical description of the isotropiza- 
tion of the particles momenta as the system evolves. Dif¬ 
ferently to the previous works we consider the balance 
between pressure anisotropy sources and isotropization 
effects, which simplifies dramatically the solution of the 
problem and results in a problem that can be addressed 
both numerically and analytically. Nature is obviously 
far more complex than a constant-A model, but still 
any approximated model that can be studied by means 
of both analytical and numerical methods may be ex¬ 
tremely usefulon the understanding of the basic dynam¬ 
ical properties related to the dynamics of collisionless 
plasma fluids. 

6.1. Magnetic field characteristic lengthscales 

Even if the amplification of the magnetic field, from 
seed fields up to equipartition levels, could be ex¬ 


plained by the collisionless plasma instabilities, one must 
still consider its characteristic lenghtscale. The zero¬ 
dimensional anaytical models cannot account for this 
properly. Still, it is reasonable to consider the typical 
lengthscale to be that of maximum growth rate. Indeed, 
from the numerical simulations presented in this work, 
it has been shown that the turbulent flows cannot mod¬ 
ify the magnetic field structure in timescales at which 
the amplification occurs. The spectral distribution of 
the magnetic field obtained from the simulations peak at 
the smallest scales possible, those at which the numeri¬ 
cal dissipation overtakes the field growth. The finer the 
resolution the smaller is the characteristic lengthscale. 

In the intergalactic medium this would be given by the 
natural process that can overtake the pressure anisotropy 
instabilities, particle collisions. With a typical mean free 
path of ^ few kpc, one should expect the magnetic field 
energy to be concentrated in fluctuations with approxi¬ 
mately this lengthscale. However, magnetic field corre¬ 
lation lengths as l arge as ^Mpc has been inferred from 
observations (e.g. iFerreti et al.ll201^ . posing the chal¬ 
lenge not only on how to amplify the magnetic field but 
also on how to distribute the energy from the small to the 
large scales. In our firehose unstable CGL-MHD models, 
we obtained power spectra which are steeper than the ex¬ 
pected for the kinetic phase of the small-scale dynamo. 
This demonstrates that the instabilities are indeed the 
dominant process in place in these models. On the other 
hand, depending on the level of pressure anisotropy it 
also results in weaker large scale fields, if compared to 
the SSD estimates. Given that most of the energy of the 
magnetic fluctuations are placed at the smallest scale in 
the gyrotropic fluid approximation (i.e. Amfp), one could 
estimate (see Section 3) the magnetic field intensity at 

large scales as Bl ~ (|A|p_l)^^^ (AmfpT”^)- For n ~ 
10“^cm“^, kT ~ leV, one finds Hmpc ~ V^IO'^^G. 

The answer to this problem may be standing in the 
properties of the turbulence itself. As discussed earlier 
in this paper, turbulence stretches and folds the magnetic 
held in a process that drives the turbulent dynamo. This 
process is dominant at small scales, and is quenched in a 
given scale when the magnetic energy reaches equiparti¬ 
tion with the turbulent one. In the collisionless plasma 
the instabilities are responsible for the growth of mag¬ 
netic held at the small scales. Even though turbulence 
cannot stretch the held lines at the smallest scales it 
could, in principle, do it at the largest ones depending 
on its strength. This is clearly shown in FigslH El and 
0 where the evolution of the turbulent Alfvenic Mach 
number is always larger than unity, even after the sat¬ 
uration of the magnetic held. This because the veloc¬ 
ity dispersion is dominanted by the large scale eddies, 
where equipartition was not reached. For these particu¬ 
lar models, the turbulence is responsible for transfering 
the magnetic held energy from small to large scales. This 
process should be even more important if the pressure 
anisotropy is quenched (this is not the case of the simu¬ 
lations, in which we have kept it constant). The quench¬ 
ing of the pressure anisotropy, followed by the transfer 
of magnetic energy from small to large scales by super- 
Alfvenic turbulence could therefore explain the observa¬ 
tions. In other words, as shown in the numerical simula¬ 
tions, the GGL-MHD instabilities are presumably domi- 
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nant over the kinematic phase of the small-scale dynamo, 
but possibly not for the later nonlinear phase (once the 
instabilities are quenched). 

Another possible solution for this problem would be 
related to the formation of large structures of matter 
(e.g. massive galaxies, groups and clusters). It is quite 
clear now that the large scale magnetic field is not pri¬ 
mordial, but has const antly evolved with the dynamical 
evolution of U niverse (iJedamzik &: Siglll2011ll . This was 
also pointed bv ISchlickeised l 20051 cf. references therein), 
in the context of small-scale fluctuations of Weibel in- 
stability. The col lapse of gas into galaxies or clusters 
(|Vazza et al.ll201^ . as well as large-scale shears, could 
in principle amplify coherent components of the mag¬ 
netic field. If correct, in such scenario magnetic energy 
would have been provided dominantly at small scales, 
in the post-recombination era, but reasembled into large 
scale coherent structures (> lOOkpc) following the re¬ 
cent structure formation of the dark and baryonic matter 
(z < 10). 


• for very weak initial seed fields, the early evolution 
of the magnetic field is dominated by collisions, and 
the turbulent dynamo should be responsible for the 
initial amplification of the field; 

• once the field grows and fli > i/n, a transition from 
an isotropic to a gyrotropic plasma occurs and in¬ 
stabilities would then be driven. A fraction of the 
“free-energy” available (|A|p) would then be trans- 
fered to the magnetic field as the system evolves, 
up to saturation; 

• saturation occurs at /3 ~ 2/|A| value, which is 

reached in a timescale Tgr oc with ijj = 

1/2 or 3/2, for |A| 1 and <C 1, respectively. 

This analytical model was then compared to a number 
of CGL-MHD numerical simulations, with different ini¬ 
tial and turbulent conditions, from which the evolution 
of the magnetic field from the seed field, up to saturation, 
was followed. The main results obtained are as follows: 


7. SUMMARY 

In this work we studied the magnetic field amplifi¬ 
cation process in turbulent collisionless plasmas. In 
such plasmas, parallel and perpendicular pressures (with 
respect to the local magnetic field) may become so 
anisotropic that instabilities may become dynamically 
important. It has been pointed in previous analytical 
models that such instabilities, acting together with a tur¬ 
bulent background, may accelerate the amplification of 
the magnetic field up to near equipartition levels. Such 
models would then have dramatic implications on our 
understanding of the magnetization of the early Uni¬ 
verse (post-recombination), given that the intergalactic 
medium then becomes collisionless. These models pro¬ 
vide timescales and saturation estimates for the magnetic 
energy density, but have never been confronted by full 
three-dimensional numerical simulations. The aims of 
this work were to revisit the analytical zero-dimensional 
models of the magnetic field evolution in turbulent col¬ 
lisionless plasmas and to provide a number of full three- 
dimensional numerical simulations that could be directly 
compared to these analytical estimates. 

As long as pressure anisotropy is kept constant during 
the evolution of the system a novel analytical solution 
is provided for the magnetic field evolution in systems 
subject to firehose instability. As discussed before, such 
assumption is not too far from reality, given that different 
processes act together to increase and decrease pressure 
anisotropy, and some sort of equilibrium value for A, in 
unstable regime, may exist (specially for weakly magne¬ 
tized plasmas). We find that: 


• the explosive growth of magnetic field energy den¬ 
sity is observed in the firehose unstable models, be¬ 
ing the amplification timescales approximately pro¬ 
portional to lAj*^, as was predicted in the simple 
analytical model; 

• there is no clear dependency between Tgr and Bq in 
the simulations, in contradiction to the linear anti¬ 
correlation derived in the analytical model. This is 
explained in terms of Bq being the field intensity 
at n ^ Vii transition, rather than the initial seed 
value. Also, because the code is unable of evolving 
the Larmor frequency as a function of B; 

• the power spectrum of the magnetic field peaks at 
the smalles scales not dominated by the numerical 
diffusion, and the transfer of energy to large scales 
must occur aftewards, by turbulence itself, once 
saturation occurs. 
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